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Abstract 

We study the impact of flavours on the efficiency factors and give analytical and numerical 
results of the baryon asymmetry taking into account the different charged lepton Yukawa con- 
tributions and the complete (diagonal and off-diagonal) L to B — L conversion A matrix. With 
■ this treatment we update the lower bound on the lightest right-handed neutrino mass. 
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^ ; 1 Introduction 

I The minimal leptogenesis scenario [T] is based on the seesaw (type I) mechanism, consisting of the Stan- 

dard Model (SM) extended by 2 or 3 right-handed (RH) Majorana neutrinos with hierarchical masses. The 
lightest RH neutrino, produced by thermal scattering after inflation, decays through out-of-equilibrium 
processes that violate lepton number, C and CP symmetries. These processes induce a dynamical produc- 
tion of a lepton asymmetry, which can be converted into a baryon asymmetry through {B + i)-violating 
sphaleron interactions. In this context, several studies investigated this possible explanation of the 

C*~) \ baryon asymmetry of the Universe (BAU) and the different analyses have led to constraints on neutrino 

physics from the requirement of a successfull leptogenesis. For instance, a lower bound on the reheating 
temperature Trh (see e.g. [H |4j |6]) and an upper bound on the absolute scale of light neutrino masses 
(for example in [7]) have been derived. These studies have been performed using the so called "one-flavour 
state" approximation. It is becoming well known that the explanation of the BAU from a successfull ther- 
mal leptogenesis has to be revisited when the mass of the decaying right-handed neutrino that produces 
the lepton asymmetry is M^ 1 < I0 12 GeV [SHU]. In this case, the Yukawa couplings of the charged leptons 
affect the dynamics of the Boltzmann equations (BE) and one cannot use the usual one-state-dominance 
approximation (one flavour), but the "flavoured" dynamical equations in the derivation of the BAU must 
be taken into account [5]-|I4j. 
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2 Flavours in leptogenesis 

We consider the SM Lagrangian with three heavy Majorana right-handed neutrinos iVj, i — 1,2,3. The 
RH neutrinos couple to the left-handed (LH) ones through the complex Yukawa coupling matrix, A. The 
small neutrino masses generated by the seesaw (type I) mechanism are given by : 

m v = v 2 U T X T M~ l XU , (1) 

where U is the PMNS mixing matrix, v is the vacuum expectation value (vev) of the Higgs field (v ~ 174 
GeV), and M. is the 3x3 diagonal Majorana mass matrix. We assume a hierarchical spectrum for the right 
handed neutrino, Mn 3 ^> Mn 2 ^> Mn 1 , and consider that the lepton asymmetry is produced by the decay 
of the lightest RH neutrino N\. Then, in the context of leptogenesis, the baryon asymmetry is obtained 
by partial conversion of the leptonic asymmetry via sphaleron interactions. With the correct treatment of 
flavoured leptogenesis, this conversion reads [9] : 

Y B ^§J2 Y ^- ( 2 ) 



Yaq is the B/3 — L a asymmetry in the lepton flavour a, which is conserved by sphaleron interactions and 
transmitted to a baryonic asymmetry 1%. 

Defining the variable z = M%/T, the BE governing the abundance of RH neutrinos Y^ 1 , and the asymmetry 
Y Aa are dJ : 

Y^(z) = -k(D(z) + S(z)){Y Ni (z)-Y*1(z)) , (3) 

YL(z) = -e aK (D(z) + S(z))(Y Nl (z)-Y^(z))-K a W(z)Y / A' afi Y A p(z) , (4) 



where Y^ q is the thermal population of the lightest RH neutrino Ni 



W^f^I*^' (5) 
with g± — 106.75 in the SM. The CP-asymmetry generated by iVi in the flavour a is given by [ID] : 

where <? is the usual loop function |15j . 

In eqs. (|3I4|) , the wash-out parameters have been factorized out, and are defined as follows: 

- AlaA lQ = — , (,g 
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with m* the equilibrium neutrino mass, m* ~ 1.08 x 10 _3 eV. These parameters exhibit the out-of- 
equilibrium condition on the decay of the right-handed neutrino: the decay process is out-of-equilibrium 
when the decay rate is slower than the Hubble expansion rate at the temperature Mm x '- V < H{Mn x ) and 
thus k < 1. The regime where k <^ 1 is called the weak wash-out regime, in which case the inverse reactions 
involving the thermal scatters are rather slow and do not efficiently wash-out the lepton asymmetry. 
On the contrary, k 1 is the strong wash-out regime, where the inverse reactions strongly wash-out 
the asymmetry. Depending on the value of this wash-out parameter, analytical approximations of the 
production and wash-out terms allow us to derive semi-analytical formulae for the baryon asymmetry, as 
will be see in the next section. The processes we take into account in eqs. (|3l4p are decays and inverse 
decays labeled D(z), and AL = 1 scattering, S(z). Notice that in eq. ((H), CP violation in scattering is 
also taken into account, and thus S(z) further contributes to the production of the Ya q asymmetry. In 
this study we neglect AL = 2 scatterings (except the real intermediate states already substracted) that 
are negligible as Af^/lO 14 GeV CO.Ixk [TT], as well as well as scatterings involving gauge bosons. The 
thermally averaged decay rate is given by: 

where K n (z) are the modified Bessel functions of the 2 nd kind. The Higgs- mediated scatterings in the s— 
and t— channel contribute to the production of the asymmetry, as well as to the wash-out term. Their 
effects can be parametrized by two functions fi(z) and f2{z): 
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(10) 
(11) 



where a/j(z) ~ log ( j^Tf\ ) — l°g(jrj) an d — 1 (2/3) in the case of weak (strong) wash-out regime. 
The wash-out term W(z) contains a part from the inverse decay and a part from scatterings [S] and is 
given by: 

W(z) = W ID (z)f 2 (z) , W ID (z) = \z z K x {z) . (12) 
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The function fi(z) (/a^)) parametrizes the effect of the scatterings in the production (wash-out) factor, 
and the r.h.s. of eqs. (|10lll[) comes from high-temperature approximations of the reduced-cross sections, 
when the scattering effects are fully relevant. In the low temperature regime, scatterings become negligible, 
so the functions /i and /2 tend to unity. The matrix A a p Q| depends on which charged lepton interactions 
are in thermal equilibrium, and parametrizes the conversion of the leptonic asymmetry into a B/3 — L a 
asymmetry according to Y a = A a pY A p. If the temperature at which leptogenesis occurs, Mi, is below 
10 9 GeV, interactions involving charged /-t and r couplings are fast compared to the Hubble expansion rate, 
and are therefore in equilibrium. Thus, fi and r flavours have to be treated separately and so the electron 
flavour is also distinguishable. Then one has [TT] : 

/ -151/179 20/179 20/179 \ 
A = -A' = 25/358 -344/537 14/537 . (13) 
\ 25/358 14/537 -344/537 / 

For M\ between 10 9 GeV and 10 12 GeV, only the charged t Yukawa interactions are in equilibrium. The 
interactions involving the e and /i flavours are slower than the expansion rate, so that those flavours are 
indistinguishable, and the decay of N\ will generate a Y e+ll asymmetry. The "flavoured" asymmetries are 
then reduced to (Y T -Y e+fJ- ), and the B — L <-> L conversion matrix reads : 



_ ( -417/589 120/589 
1 30/589 -390/589 



(14) 



In a recent work [T3] , it has been argued that the interaction rates involving the charged Yukawa couplings 
should be fast compared to the interactions involving the decaying N\ in order to have sufficient time to 
project the produced lepton asymmetry onto flavour-space. It has been derived that M^ should be below 
5 x 10 11 GeV for the tau- Yukawa to be in equilibrium, hence projecting the lepton asymmetry on the 
(Y T -Y e+fJi ) space. This point will be discussed in section 4. 
A formal solution of eq. ((3]) for the B/3 — L a asymmetry is given by: 



Y Aa (z) = -e Q K / dxD(x)f 1 (x)AN 1 (x)e~ K " A ^^ dyW ^ (15) 
- K a f dxW(x)Y AfJ (x)e~^ A '™ i: dyW ^ , 

where ANi(z) — (Yjvi(z) — Y^(z)~) is the departure from thermal equilibrium. The first term in eq. (fT5|) 
had been estimated for a vanishing initial N\ abundance, Ni(zi„) — 0, and for an N\ abundance initially at 
thermal equilibrium Ni{z in ) = Ni(zi n ) eq , in Refs. 0[TTJ[T2]. The second term, which has been neglected 
in these previous studies, is responsible for the interdependency of the flavours through the off-diagonal 
matrix elements A' a g. This term drives a new contribution to the B/3 — L a asymmetry and can be relevant. 
In some cases the latter is in fact the dominant contribution, as we will later see. 

We parametrize, as usual, the final asymmetry Y Aa in terms of efficiency factors that contain all the 
dependency on the wash-out factors k, n a . Those cfficncics r/ a arc defined by: 

Y Aa = -e aVa Y^{T-»M Nl ) 

~ -3.9xl0- 3 £ Q (^+< d ) . (16) 

The first term rj d has been derived in jTTJ [TJ] , and its expression will be presented in the next section. The 
second term, ?7™ d , arises from the non-diagonal conversion of a leptonic flavour, say Lf3 1 into the B/3 — L a 
direction and its effect is also studied in the next section. It is clear from eq. (|15l) that the efficiency r) a 
of the process depends on the individual wash-out parameter k q , but weighted by the factor A' aa arising 
from the B — L <-> L conversion. We then define k a = A' aa n a (and consequently k — ^ Q k a ) as the "real" 
wash-out parameter, and thus r\ a = r](R a ). 
The baryon asymmetry reads: 

Yb = P Fa « - - 1 ' 26 X 10_3 J2 e « Vc* - With Va=vt + Vl d ■ (17) 



For convenience we use A 1 = —A in the BE so that the diagonal elements are positive. 
2 We will refer to the case of a vanishing initial abundance as a dynamical case, as the population of N\ is created 
dynamically by thermal processes. The case of a iVi initially at thermal equilibrium will be refered as thermal case, even if 
it recquires a non-thermal production mechanism. 



3 



The baryon asymmetry is, as the lepton asymmetry, the sum of the diagonal term proportional to A aa ~ 1 
and of the off-diagonal term proportional to A a p_ p^ a ~ 1/10. The latter contribution will be negligible 
for the baryon asymmetry, but will strongly modify the individual lepton asymmetries. 

3 Efficiency factors 

Here we proceed to numerically solve the BE (eqs. ()3l4j) ) for different configurations of the individual CP 
asymmetries and wash-out factors, for distinct flavour "alignments" . 

3.1 Study of the efficiency rf 

In this part, we neglect the off-diagonal part in the last term of eq. and solve the BE (eqs. (13l4j) ). 
Depending on the initial conditions and on the strenght of the wash-out, several analytical approximations 
can be derived for -rf (eq. (fT6|) ). 

3.1.1 Vanishing initial Ni abundance 

In the case where the population of Ni is dynamically generated, i.e. Ni(zi n ) — 0, one can derive the 
expression for rj d in different wash-out regimes [llj : 

• all flavours in the strong wash-out case: n a ^> 1 

/ 1 \ 116 

^- 3 - 5 fe) ' (18) 

• all flavours in the weak wash-out case: k„ « 1 , 

if{k a ) ~ 1.4 /c a ii . (19) 

However, within the choosen range of wash-out parameter, we find that the efficiency is better fitted 
with ri d {k a ) ~ 0.4 k a \fk , as can be seen in figffl 

• some flavours in strong n a 3> 1 and some others in weak np <C 1 wash-out regimes. In this case the 
efficiency for the flavour a is given by eq. 1|18|) and the efficiency for the flavour j3 is given by: 

r, d (kp) ^ 0.3 hp . (20) 

The efficiency rj d is then obtained by simple interpolation between these three generic cases. For the sake 
of illustration, we choose 3 representative cases: 

• Case a): all the wash-out parameters K a are equal. 

• Case b): some flavours (e.g. 0) are weakly washed-out with Kp — 5 x 10~ 2 . 

• Case c): some flavours (f3) are stronlgy washed-out with up = 30. 

We checked the validity of those expressions in the considered range of wash-out parameters, n a between 
10~ 2 and 10 2 . It is interesting to notice the dependence of r) d on the total wash-out parameter k in 
eq. (fT9"|) . From this term, a flavour a that is weakly washed-out will be sensitive to the wash-out of the 
other flavours implying that there is a correlation of the flavours. This can be seen in figure [lj where we 
represent the efficiency of a given flavour (a) as a function of the respective wash-out parameter n a , for the 
three representative cases a), b) and c) discussed above. The democratic scenario, case a), which is similar 
to the one-flavour approximation, is represented in red and the misaligned cases b) and c) are represented 
in blue and green, respectively. For comparison, we also represent the analitycal estimates (dashed lines) of 
the efficiencies, eqs. (|18M20p . We see that the agreement between the numerical and the analytical results 
is very good. 
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Figure 1: Efficiency r] d (k a ) in the dynamical case, for specific values of Kp, f3 ^ a. Case a) (red 
curves): up = K a ; case b) (blue curves): np — 5 x 10 -2 ; case c) (green curves): up = 30. In each case, the 
solid lines represent the numerical computation and the dashed ones the results of the analytical approximation. 
In the left panel, the upper (dotted) curve corresponds to the analytical estimates of the efficiency in the weak 
wash-out regime eq. (|19|l . 

Firstly, in the case where the right handed neutrino population is created by inverse decays and scat- 
terings, the wash-out factor which parametrizes the strenght of those thermal production, has to be non- 
negligible in order to produce a sufficient amount of iVi. Therefore, the efficiency is maximized in this 
dynamical case for n a ~ 1 with r\™y* — 0.2. 

Secondly, and this is the main point here, we see that the efficiency of the process when the flavour 
a is weakly washed-out does indeed depend on the strenght of the wash-out of the other flavours. For 
example, for n a ~ 5 x 1CP 2 (and assuming A' aa ~ A'gg ~ 1), we roughly have 1$ ~ rja ~ 3 x 1CP 3 and 
r$ ~ 1.5 x 1CT 2 . 

The enhancement 77°' comes from the fact that n ~ 10 but still K a < 1, and the eq. ([2"0|) applies. For 
cases a) and 6), eq. ()19p applies, and we have an extra supression from the factor k ~ 10 _1 . Looking at 
the strong wash-out regime, K a > 1 , we see that the efficiency of the flavour a does not depend on the 
wash-out of the other flavours. 

3.1.2 Thermal initial abundance 



In the case where the population of N± is initially in thermal equilibrium, Ni(zi n ) — Ni(zi n ) eq , the 
computation of the efficiencies is modified. Following [51 [H], one has for the efficiency factors: 




icr 2 icr 1 1 10 10 2 itr 2 icr 1 1 10 io 2 io -2 io _1 1 10 10 2 



Figure 2: Efficiency rj d {k a ) in the case of a thermal initial N\ abundance, for specific values of Kp, f3 ^ a. Line 
and color code as in figure [T] 

The study of the thermal case is shown in figure where we again consider the 3 cases a), b) and c), 
with the same colour code as in figure Q] The difference from the case of a vanishing initial Ni is striking 
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in the weak wash-out regime: besides the obvious behaviour of rj that is maximized for small wash-out, 
^ther — 1 for k q <C 1, we see that the effect of flavour is negligible, except for K a — 1, where a small 
distinction of the cases a), b) and c) is possible. In the strong wash-out regime, the individual efficiencies 
do not depend on the wash-out of other flavours, as in the dynamical case. Furthermore, in this strong 
wash-out regime, there is no distinction between the thermal and the dynamical cases, as can be seen in 
the left panel of figure 2] which is included in the discussion of the next section. We also see that the 
agreement between the analytical result for the efficiency given in eq. (|2ip and the numerical result is very 
good. 
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3.2 Study of 77 

Now we study the effect of the non-diagonal terms A a p ,a^f3, that exhibit the interplay between different 
flavours. This interplay can be seen in the last term of eq. (fT5|) : 



- Y2* = 3.9 x Itr 3 e aV f = n a J2 f dx W(x)Y A0 (x) e -"«^» . (23) 

One can find an approximate expression for this term: considering that the variations of Y A „ (as) are 
negligible compared to the rest of the integrand, one can approximate it to its final value Ya 3 (as) — Ya (00), 
and can thus be factorized out from the integral. Another approximation is to consider that Yao is mainly 
generated by the diagonal part, that is we neglect the effects of 0{A^ ld ) that are corrections of the order 
of a few percent. We obtain: 

3.9 x 10- 3 e a r?2 d * * a ]T A' a0 Y^ f c (k a ) , (24) 

where k a = n a A' aa and f c (k a ) is given by: 

f c (k a ) = I dx W(x) e -~ K «f~ dvW M (25) 
1 

1 + 0.8x kl 



L3 , , n o gl .X7 • (26) 



The total efficiency of a given flavour is the sum of the contribution from the diagonal part of A which, 
contains slight contamination from the other flavours, and from the non-diagonal part, that will be re- 
sponsible, as we will see, for a huge modification of the total efficiency, in some cases becoming dominant 
compared to the diagonal contribution, 

j](k a ) = i] d (k a ) + n a f c {k a ) A a/3 — rf{Kp) ■ (27) 

The effect of the non-diagonal part on the asymmetry Y& a depends on the wash-out and on the CP 
asymmetries of the different flavours. For example, if we consider the flavour a, the asymmetry produced 
by the diagonal part proportional to r] d depends on the strenght k a . If k a 3> I, then Ya q is strongly 
washed-out and therefore is too small, as can be seen in figure [T] for k q > 100 we have r] d (k a ) < 10~ 3 . 
Now consider the non-diagonal part: it depends on the strenght of kp. If the wash-out of the flavour (3 is 
weak (or even mild), then rj d {kp) will be close to its maximal value, and then one has in this case: 

V (k a ) cc ^-4 max x (lO^-lO- 2 ) , 

the value of ?7 max depending on the inital conditions. Typically, this effect can, in favourable situations, 
drive up the asymmetry Y\ a by one or two orders of magnitude as can be seen in the different panels of 
figure [3] There we represent the efficiency factor of a given asymmetry r\ a as a function of the wash-out 
parameter K a , for the different alignments of the flavours. We clearly see that the off-diagonal terms of 
A modify the efficiency rj a in the strong- wash-out case n a 3> 1. Another consequence of this dependence 
of ?7™ d on the wash-out of the other flavours is illustrated in figure [H We see that when we neglect the 
off-diagonal terms there is no distinction in the strong wash-out regime between the dynamical and the 
thermal case. That is, in the strong wash-out regime, the efficiency is independent of the thermal history 
of the decaying Ni. However, when we include the off-diagonal terms, we see that for the case c), the 
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Figure 3: 

Influence of the other flavours on the efficiency r/(ii a )- Upper panels: case of a vanishing initial abundance. Lower 
panels: case of a thermal initial abundance. The curves represent different inputs for the wash-out parameters Kp, 
(3 7^ a: for the case a) (red curves) we have a democratic scenario where Kp = K a , while for the case b) (blue curves) 
Kp = 5 x 10 -2 and for the case c) (green curves) Kp = 30. The solid lines represent the numerical computation and 
the dashed ones the results of the analytical approximation. 
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Figure 4: 

Influence of the other flavours on the efficiency rj{k a ). The left panel represents the diagonal contribution to the 
efficiency, whereas the right panel represents total, the diagonal and non-diagonal contributions. We represent the 
case of a zero and thermal initial Ni abundance, for the three cases a), b), c). Line and color code as in figure [T] 

two efficiencies do not coincide anymore in the strong wash-out regime, as each f](k a ) is related to 77^(^,3), 
the latter strongly depending on whether the Ni was initially at thermal equilibrium or had a vanishing 
abundance. This clearly illustrates the effect of flavours in leptogenesis. 

The above discussion is based on a strong approximation (Y/^ /3 (x) ~ Ya^(oo)), which was made in 
order to quantify the off-diagonal effect in eq. (|23[) . However this excessively naive approximation does 
not describe all wash-out and CP asymmetry configurations. Furthermore, the dynamics of the flavours 
(3 7^ a is neglected and in particular, in the case of a strongly washed-out flavour /3, this approximation 
cannot be used. We thus solve the BE, eqs. (J3|), {3|, including the off-diagonal terms of the A matrix and 
illustrate in figure [5] contours of the ratio Y e+f ^/Y T (absolute value) in logarithmic scales, for the thermal 
and dynamical cases, as function of the ratio of the flavoured CP asymmetries e e +^/ e T and of the flavoured 
wash-out parameters K e+M /K r . At first sight, we see that Y e+il > (<)Y T for e e+M > (<)e T . 

Looking in more detail, we see that the wash-out parameters influence the former statement. For 
example, in the thermal case, for k £+ ^/k t ~ 0.1 and e e + M /e T — 0.5, we have Y e+fl ~ 2Y T . 
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Figure 5: Contour plots of the ratio of the individual flavoured asymmetries, plotted as a function of the ratio 
of wash-out parameter versus the ratio of CP asymmetries. We place ourselves in a 2 flavour case with Mni = 
5 x 10 9 GeV, € e +fj, + £t = 10~ 6 and re e +/j + k t — 10. The left (right) panel stands for the dynamical (thermal) case. 
On both graphs, the black area denotes the equality Ye+y, = Y T . Moving to the left of this countour, Y e +^ < Y T 
whereas to the right we have ie+p > Y T . 

3.3 The baryon asymmetry 

Summing-up the contribution from the diagonal and off-diagonal parts of the A matrix, we finally obtain 
the baryon asymmetry: 

Y B = ^^r Act ~-1.26xlO- 3 ^e ct ? 7a . (28) 

a a 

For e = ^ Q e a ^ 0, one can define an efficiency for the baryon asymmetry t\b such that: 

Y B = -1.26 x 10- 3 er) B . (29) 

The baryon asymmetry will be the sum of the individual lepton asymmetries. Therefore, the baryon 
asymmetry, and hence the efficiency r/s, will depend on the alignment of the flavours. We illustrate this 
point in figure [51 where the efficiency t\b is plotted as a function of the wash-out parameter n a , for cases 
a), b) and c). The solid lines represent the diagonal contributions, and the dashed ones represent both 
diagonal and non-diagonal contributions. As the baryonic efficiency can be defined only for e ^ and 
is strongly dependent on the flavoured CP violation e a , we set e a = 8 x 10~ 7 for all flavours. First, we 
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Figure 6: Effects of the off-diagonal terms of the A matrix on the total baryon efficiency r\s- 
consider only the diagonal contributions. As expected, we see that the efficiency strongly depends on the 
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flavour alignment. The democratic scenario (case a), in red) shows the one- flavour approximation-like 
behaviour. The flavours play a full role in the misaligned case, and especially in case b) (in blue), where 
we have Kp = 5 x 10~ 2 . Then the efficiency for the flavour(s) j3 is (are) close to its (their) maximum, 
depending on the thermal history of Ni. If the flavour a is strongly washed-out, the baryon asymmetry is 
mainly composed of the flavour(s) /3, and is thus weakly washed-out, even if the total wash-out is strong. 
This is the very effect of taking into account the flavours in leptogenesis. We also see that the baryonic 
efficiency, as well as the leptonic one, depend on the initial abundance of Ni in the strong wash-out regime. 
Now, let us consider the dashed lines, which represent the sum of the diagonal and off-diagonal terms of A. 
We see that the off-diagonal terms account for percent corrections. A more detailed analysis of their effect 
on Yb is shown in figEl where we represent contours of Y B tot ^ /Y g laeon in the interesting two flavour-case. 
We choose the individual CP asymmetries to be equal, e e+Al = e T (= 1CP 6 but the actual value is not 
relevant). We see that the off-diagonal terms affect the baryon asymmetry up to 40% in the dynamical 
case in both ways (increasing and decreasing). On the other hand, the thermal case is only enhanced, up 
to the same amount. Effects of the off-diagonal terms are particularly important in the strong wash-out 
regime. In some particular cases the non-diagonal terms also have a relevant role, namely in the democratic 




Figure 7: 

Effects of the off-diagonal terms on the baryon asymmetry: contour plot of y^ otal /Yg iag ' for fixed CP asymmetries 
6e+/4 = £t and varying wash-out parameters. The left (right) panel stands for the case of a vanishing (thermal) 
initial Ni abundance. 

scenario, where k a are equal, and where the total CP asymmetry is zero ^ a e a = 0. Then ignoring the 
O(A) effect will lead to a vanishing baryonic asymmetry, as ^ Q e a rfe = [17] • On the contrary, taking 
into account those effects avoids the latter cancellation. For example, in the case of three distinguishable 
flavours with the specific alignment in which k a — kp — kg and e a + ep + eg = e Q (l + b + d) — 0, one finds 
for the baryon asymmetry : 

Y B = -1.26 x l(r 3 e Q T] d (k a )\ K ai fc(k ai )A ai a + b ^2 k Pi f c (kp 1 )Ap 1 p + k Si f c (kg 1 )Ag 1 g \ ■ 

In order to maximize the above function, we take a = fi, f3 — e and 5 — t for a positive value of b (in the case 
of a negative value of 6, a = fx, (3 = r and S = e). We then get ~ e^^k^) f^k^n^ x4.5 x 10~ 5 (1 + 2b ), 
which is non-zero for b ^ —1/2. 



4 Lower bound on 

When the flavours are taken into account in leptogenesis, the modification of the asymmetry, combined 
with the change in the efficiency factor may have an impact on the lower bound of the mass of N±. From 
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the bound on each individual CP asymmetry [10] , 



< 



16ttv 2 



(30) 



one has 



\Y B \ 1-26 x 10- 3 J2 



< 1.26 x 10 



from which a lower bound on Mm, is derived 



_ 3 3 M Nl m max 
16ttv 2 



M 



JVi 



> 



16tt 



M 



iVi 



> 7.1xl0 8 GeV 




(31) 

(32) 
(33) 



Since the lower bound on Mjsf 1 is inversely proportional to the efficiency r\ a , it will therefore depend on the 
thermal history of the decaying right-handed neutrino. In the case where N\ are produced by scatterings, 
the efficiency is maximized for a wash-out K a ~ 1, where rj a ~ 0.2. In the case where N± are non-thermally 
produced, the efficiency is maximized to 1 for a very weak wash-out k„ « 1. The lower bound will depend 
on the alignement of flavours, and in the democratic case one has: 



M Nl 



> 



4.1 x 10 s GeV in the thermal case 
2.5 x 10 9 GeV in the dynamical case . 



(34) 



This bound is close to the one derived in the "one flavour approximation", where M^ 1 > 2.1 x 10 s GeV 
in the thermal case and Mn x > 1.05 x 10 9 GeV in the dynamical one [7J. Besides flavour effects, the 
difference between the lower bounds of the flavoured and unflavoured cases comes from a different factor 
in the B — L <-> B conversion. Indeed, in the one flavour dominance, the Davidson-Ibarra bound reads [6]: 



e < 



3 M^rrir, 
8^ 



(35) 



and the conversion from sphalerons is Yb = 28/51 Yi,. Therefore the lower bound on Afjy-j in the unflavoured 
case is ~ 1/2 x 12/37 x 51/28 ~ 0.3 times the flavoured one. 

An implication of this bound resides in the well-known conflict between the reheating temperature and 
leptogenesis. Indeed, Trh should be above M^ 1 in order to avoid the erasing of the lepton asymmetry. 
In view of the above estimates, Trh ^ 4 x 10 9 ' 8 ' GeV in the dynamical (thermal) case. In this case, the 
inclusion of flavours does not really help. In the regime of strong wash-out, K ^> 1, where the effective 
neutrino mass is close to the mass inferred from atmospheric oscillations, the situation changes. In the 
one-flavour approximation, the efficiency t](k) oc k , therefore M^" 11 oc k, and increases with the wash- 
out, and so does the reheating temperature. In this strong wash-out regime, we roughly estimate Trjj > 
M/VjVlO ~ 9 x 10 8 GeVK, k> 1 [|. Including flavours, one generically has non-alignment of the individual 
asymmetries, and one can have a strong total wash-out, with some weakly washed-out flavours. If for 
example (see fig. [SJ, k ~ m atm /m* ~ 45 with Kp ~ 40 but n a ~ 5, one has M^ n ~ 2 x 10 10 GeV for 
the flavoured case and 8 x 10 10 GeV for the unflavoured case. Hence, including flavours, the reheating 
temperature in this strong wash-out regime is lowered by a factor ~ 4, and one roughly has Trjj > 
2 x 10 9 GeV. 



4.1 Numerical results 

We numerically solve the set of coupled BE (eas.(|3ttp) and obtain the allowed parameter space. The input 
parameters of the BE are the Yukawa couplings A, which define the wash-out parameters and the CP 
asymmetries. These have been built using the Casas-Ibarra parametrization |16j . In this parametrization 
the matrix A reads: 

A = Mj/ 2 RmWrfv- 1 , (36) 



10 



U being the PMNS matrix parametrized by three angles and three phases. The solar and atmospheric 
angles are well measured, whereas the #13 angle (Chooz angle) is only upper-constrained. The three CP 
violating phases are yet undetermined. The matrix R is a 3 x 3 orthogonal matrix depending on three 
complex angles. We impose a normal hierarchical spectrum, both for the light and for the heavy neutrino 
sectors, so that the only free parameters that enter in the mass matrix are the lightest neutrino masses mi 
and Mn x ■ Therefore, we have 12 independent variables that are not experimentally determined. 

In our numerical computation, we scan over the whole parameter space by randomly choosing the 
free parameters. Moreover, we impose a perturbative limit A33 < yt (yt being the top Yukawa coupling) 
that will upper-constraint Mjv i: leading to the upper bound < 4 x 10 11 GeV. From [14] and in 

order for the flavour to be relevant in leptogenesis, an upper bound on MjVj can be derived, namely 
Mjvj < 5-8 x 10 11 GeV. Here, we only consider cases were the interactions involving the charged lepton 
Yukawas are faster than the inverse decay, so that flavours are fully relevant. Finally, by imposing the 
obtained baryon asymmetry to be in the experimental ranged, 5.2 x 10~ 10 < 7.04 x Yg < 7.2 x 10 -10 , we 
represent in figure[8]the (Mn x — rh) parameter space allowed by the requirement of a sucessfull leptogenesis. 
Many remarks are in order. Firstly, the lower bound numerically derived confirms the one given in eq. (|34p . 
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Figure 8: 

Successfull leptogenesis: bound on Mi, in the case of a zero (thermal) initial Ni abundance for the up (down) 
panels. The left panels show the allowed (Mtf^m) parameter space for the one-flavour approximation, whereas the 
right panels stand for the case when lepton flavours are taken into account. The vertical lines represent y Am^ tm 
(in blue) and a/ Am^, (in green). 

and we notice that the "one flavour approximation" (left panels in figure [S]) lowers the bound compared to 
the correct flavour treatment (right panels). Secondly, comparing the dynamical and thermal cases, that 
is the up and down panels in figure [SJ we see that m can take much smaller values in the case of a thermal 
initial abundance. As this mass encodes the out-of-equilibrium condition, cf. eq. JJJ, this means that in 
the thermal case, leptogenesis occurs even for extremely small values of m, while in the dynamical case, 
this is not possible. Finally, comparing the left and right panels, one clearly sees the effect of flavour in 
the "re-opening" of the parameter space for higher values of rh. Indeed, by introducing the lepton flavour 
asymmetries, we relax the one-flavour approximation that corresponds either to a common behaviour of 
the individual asymmetries, or to the dominance of one flavour. Now, other configurations are allowed, 
and the misaligned ones widen the parameter space. The mass rh is related to the total wash-out, that is, 
to the sum of each individual wash-out. In the one-flavour approximation, a high value of rh corresponds 
to a strong wash-out and is disfavored by leptogenesis. On the contrary we observe (cf. figure O that even 



3 The experimental bounds come from cosmological constraints: the lower bound corresponds to the Big-Bang nucleosyn- 
thesis (BBN) constraint on the light species abundance, and the upper-bound from cosmic background radiation (CBR) 
constraints 1181 . The WMAP 1191 constraints are included therein. 
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if the total wash-out is strong, we can still have flavours that are weakly washed-out, hence dominating the 
baryon asymmetry and allowing a successfull leptogenesis. Thus, by the inclusion of flavour in leptogenesis 
no upper-bound on m can be derived. 

Notice that for m(mi) > atm in figure [3[9|) , the points drop below the upper-bound Mjvj — 5 x 10 11 GeV. 
Indeed, as mi > atm , m max ~ mi ~ 1112 — W3, and the upper-bound on MjVi scales as 1/mi, c.f eq. (|32p . 
In the one flavour approximation, a bound on the light neutrino mass was derived in [7], and this no 
longer holds when flavours are accounted for [10] . However, in [13j a bound on the neutrino mass scale of 
about 2 eV is derived in the flavoured leptogenesis context in the strong wash-out regime and hierarchical 
wash-out factors 1 <C K a <C Kp and equal CP-asymmetries. In this work, we impose m„ to be lighter than 
the cosmological bound ^2 m v < 1 eV and we do not explore configurations leading to higher m„. This can 
be seen in figure O which represents the allowed parameter space (M^-mi ) in different cases. The black 
points are the result when flavours are included, whereas red ones represent the one-flavour approximation. 
We clearly see that the cosmological bound is saturated when flavours are considered, and this does not 
occur in the one-flavour approximation. In figure for m\ above m a tm , the solutions have in general a 
specific flavour alignement: the flavoured CP-asymmetries are almost equal e e +n ~ £t and the individual 
wash-out factors are hierarchical 1 <C n a < 10np and the total wash-out is strong. It is well known that 
such configurations of wash-out parameters achieve sucessfull leptogenesis in the flavoured case whereas 
the unflavoured one fails. For such specific configurations, effect of the off-diagonal terms of the A-matrix 
on Yb is maximisal (c.f fig[7|) but nevertheless is only a correction without important impact. 
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Figure 9: 

Successfull leptogenesis: M\-m\ space, in the dynamical case (left panel) and in the thermal case (right panel). 
The vertical lines represent yj Am^ tm (in blue) and yl Am^ ol (in green). 



5 Conclusion 

The behaviour of individual lepton asymmetries in the case of vanishing initial N\ abundance has been 
analysed in |llj . In this study we give semi-analytical results including fine-tuning corrections that depend 
on flavour alignment. We extend the study to the case of N\ initially in thermal equilibrium, and confirm 
that in this case, when off-diagonal entries of the conversion B/3 — L a <->■ L are neglected, the efficiency 
factor for a given flavour is independent of the wash-out of other flavours. 

Independently of the thermal history of the decaying right-handed neutrino, we observed that misalignment 
of flavours can greatly enhance the baryon asymmetry, when compared to the one-flavour approximation, 
for an identical wash-out strenght.We also include off-diagonal entries of the B/3 — L a «-> L conversion 
that couple flavours among themselves. Even if this inclusion only modifies the baryon asymmetry by a 
few percent, thus allowing to safely disregard these terms for Yb computation, we nevertheless stressed 
that the individual lepton asymmetries are very sensitive to such interdependencies. Finally, we studied 
the lower bound on the N± mass and the leptogenesis allowed parameter space. We confirm the lower 
bound to be ~ 4 x lO 8 ^ 9 ) GeV for a thermal (vanishing) initial Ni abundance. We have also shown that 
the parameter space is enlarged, as the flavour (mis)alignment allows for higher values of the wash-out (or 
equivalently of m). 
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